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THE DIFFERENCE BETWEEN A DISCRETE AND CONTINUOUS 

HARMONIC MEASURE 

JIANPING JIANG AND TOM KENNEDY 


Abstract. We consider a discrete-time, continuous-state random walk with steps uni¬ 
formly distributed in a disk of radius of h. For a simply connected domain D in the 
plane, let Uh{0,-;D) be the discrete harmonic measure at 0 € D associated with this 
random walk, and w(0,-;D) be the (continuous) harmonic measure at 0. For domains 
D with analytic boundary, we prove there is a bounded continuous function <jd(z) on 
dD such that for functions g which are in C 2+a (dD) for some a > 0 


lim 

740 


J dD 9(Ouh( 0, |d£|; D) - f QD s(£M°> |d£|; D) 


I dD 


g(z)o D (z)\dz\. 


We give an explicit formula for ud in terms of the conformal map from D to the unit disc. 
The proof relies on some fine approximations of the potential kernel and Green’s function 
of the random walk by their continuous counterparts, which may be of independent 
interest. 


1. Introduction 

Let D be a simply connected domain in the complex plane with z G D. Suppose B( 
is a complex Brownian motion started at z, and Td '■= inf (t > 0 : Bf e( D} is the first 
exit time. The (continuous) harmonic measure at z is the probability measure on dD, 
u>(z, •; D ), defined by 

u(z,r-,D) = p(b z td e r), (i) 

where T C dD. See [7j for more information about continuous harmonic measure. 

We place a square lattice with mesh size h > 0 on D. Suppose S* h is a simple random 
walk in hZ 2 started at Zh where Zh is a closest lattice point to z, and To : = inf{n > 
0 : S* h ^ D} is the first exit time. Then the discrete harmonic measure at Zh is the 
probability measure on dD, u)h( z h, S D), defined by 

u h (z h ,T;D) = p(s^ e r) 

where Si^ D is a point in dD that is closest to S^ D and T C dD. It is well-known that 
Wh{ z h i S D) converges weakly to u(z, •; D) as h J, 0. We would like to study the difference 
uJh — oJdsh—^ 0. We expect this difference to be of order h, and we would like to compute 
the limit of (u>h — iS)fh as h approaches 0. 

The simple random walk is not rotationally invariant, but its limit as the lattice spacing 
goes to zero is rotationally invariant. However, there is no reason to expect the limit 
lim/qofU/). — co)/h to be rotationally invariant. Our simulations indicate that this limit 
depends on how the lattice is oriented with respect to the domain. We have not been able 
to rigorously study this limit for the simple random walk. In this paper we consider a 
random walk that takes place in the continuum and is rotationally invariant even before 
we let the lattice spacing h go to zero. For this model the ratio —ca)/h is rotationally 
invariant. 

The continuous-state random walk: The walk we study is given by S n = Sq + Ad + 
Ad + • • • + X n where Ad’s are i.i.d and Ad is uniformly distributed in the disk of radius 
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h > 0. We use P x to denote the conditional distribution of given S 0 = x ; and 

write E x for the corresponding expectation. 

Let T d — inf{n > 0 : S n ^ D}. Now the discrete harmonic measure at x is defined by 

uj h (x,T;D) = P X (S^ G T) (2) 

where St d is the orthogonal projection of St d onto dD whenever such orthogonal pro¬ 
jection is well-defined (it is if dD is analytic and h is small) and T C dD. 

Before stating our main result, we need some terminology. For z G C, we denote the 
imaginary part of z by Im(z). Let HI = (z G C : Irn(z) > 0} be the upper half-plane. For 
z G dD, we say dD is locally analytic at z if there exists a one-to-one analytic function 
/ : D := {£ G C : |£| < 1} —> C such that /(0) = z and 

/(D) flh = /(Dni). 

We say dD is analytic if dD is locally analytic at each point of dD. 

Our main result for the continuous-state random walk is the following theorem. 


Theorem 1. Suppose D C C is a simply connected and bounded Jordan domain, and 
dD is analytic. Assume 0 G D. Then there is a bounded continuous function <Jd(z) on 
dD such that 

lhl J d B9(OMUAdZ\;D)-f dD g(Ou J (0,\dZ\;D) f 
HO h J dD 

for all functions g on dD which are C 2+a with respect to arc length along the boundary 
for some a > 0. Uh is defined in ([2]), and to is defined in flT]). 


We prove the theorem in two steps. The first step is the following proposition. 


Proposition 1. Suppose D satisfies the conditions in TheoremUl Let g be a function on 
dD such that the harmonic function f on D with boundary data g is in C 2 (D). Then 


lim f dD g(IH(Q, Kl;-P) - J 9D g«MQ. W; d) 

HO h 


K 



Hd( 0 , z) 


df_ 

dn z 


(z) \dz\, 


(3) 


where Hd(0,z ) is the Poisson kernel (i.e o;(0, \dz[, D) = Hd( 0, z)\dz\) and n z is the 
inward unit normal at z. The constant K is given by 


K = + — [ (sin 2 ^ — (sin 4 #)/3 — 0 cos 9 sin 0)Ejf°l e (llm^STgd^dO, 

457T 7T J Q 

where Ejff^ 9 is the conditional expectation of {S n } n > 0 with h = 1 given S 0 = i cos 9, and 
T m := inf {k > 0 : S k i H}. 


Theorem jT] follows from Proposition |T] if we show that there is a function &d(z) on the 
boundary such that 

K f H d (0,z)-^-(z) \dz\ = f g(z)a D (z)\dz\. (4) 

Jod dn z J dD 

This follows from Proposition [3] in Section [4] and a trivial change of variables. Proposition 
[3] gives an explicit formula for ud- 


Remark 1. Suppose 'y(s), 0 < s < length(dD) is the arc length parametrization of dD. 
If 7 (s) G C 2+a and gi'y(s)) G C 2+a for some a > 0 then corollary II. f. 6 of [7] implies 
that f G C 2 (D). Note that this implies that if g satisfies the hypothesis in the theorem, 
then it satisfies the hypothesis in Proposition [7} 
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Remark 2. Proposition [7j is proved when each step of the random walk (i.e., X t ) is 
uniformly distributed in the disk of radius ofh. Actually, the same result but with different 
K holds if: the distribution of X, is rotationally invariant, X, is supported in the disk of 
radius of h, and the potential kernel for such a random walk has similar asymptotics as 
described in Lemma W Le -, a(x ) = C\ In |x| + C 2 + 0( |x| 2 ) for some constants C\ and 
C 2 ). The proof for such general X t is similar to the one of Proposition^ 

Remark 3. We believe Theorem [I] also holds for piecewise continuous functions g. 

Remark 4. The complicated expression for K is what comes out of our proof. Monte 
Carlo simulation of this expression gives K = 0.2647664 ±0.0000026. We conjecture, but 
cannot prove, that K is also given by the much simpler expression 

K = lim £^ =1 (|Im(S TH )|). (5) 

y—>00 

Monte Carlo simulations of these two formulas for K give values that agree to within 
about 10“'. 

The proof of the theorem relies on fine estimates of the potential kernel and Green’s 
function for the continuous-state random walk. We prove the scaled Green’s function of 
the continuous-state random walk is close to the usual continuous Green’s function up to 
O(h). We give an estimate up to 0(h 1+e ) for some e > 0 when the Green’s function of 
the continuous-state random walk is evaluated at points near the boundary. 

There is a close relationship between harmonic measures and Dirichlet problems. The 
continuous Dirichlet problem is to find a harmonic function / in D with prescribed 
boundary values on dD. More precisely , find / G C(D ) D C 2 (D) such that 


= U, z 


A/W = ^ + ^ = 0, 
/0) = g(z), 


G D 
z G dD. 


( 6 ) 


Of course the existence of a solution depends on the smoothness of dD and g, and 
uniqueness can be proved using the maximum principle. If D is regular and g is continuous 
then the solution of the problem ((HD can be written as 

f(z) = E[g(Bl D )\= I g(Ooj(z,\d(\;D). 

J dD 

See for example Theorem 8.5 of [T3j for a proof. 

The generator, Ah, for the continuous-state random walk is defined by 

A*/(z) = / if( z + 0 ~f( z )] d ^ ( 7 ) 

7r "' 2 B(0,h) 


where B(0,h ) denotes the disk with radius h centered at 0, and df denotes the usual 
two-dimensional Lebesgue measure. 

We divide D into two subdomains, D 2 {z G D : dist(±, dD) < h} and Di := D \ D 2 
where dist(z, dD) = inf{|^ — w\ : w G dD}. Also, we let D 3 := {z G C \D : dist(z, dD) < 
h}. The discrete problem that we consider in this paper is defined by 

^hfh (z) = 0, zeD 

f h (P) — g{P)i PtD, 

where P G dD satisfies \P — P\ = min{|z — P\ : z G dD}. Note that such P is unique 
when dD is analytic and h is small. It is easy to check that fh(z) = E z [g(ST D )] = 
J dD g(f)uh(z, |d£|; D) is a solution of (j8]), and the uniqueness follows from the maximum 
principle described in Lemma |2j 


( 8 ) 
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In the PDE literature, fh — f is usually called the “discretization error”. If one defines 
(JED using the generator for the simple random walk, then | fh — f\ — o(l) was established 
in pf] while \fh~ f\ = 0(h) was proved in pj. See pSJ, Section 23 of 0.0 and references 
therein for more about the discretization error. Since f g(f)ujh(z, |d£|; D) = fh(z ) and 
f g(f)u>(z, |d£|; D) = f(z), our theorem implies that \fh — f \ — 0(h). To our knowledge, 
all the existing results with discretization error of 0(h) assume that / G C a (D) (see [5] 
or Section 23 of 0)- Our theorem proves the same discretization error 0(h) but only 
assuming / G C 2 (D). 

The organization of the paper is as follows. In Section [2] we define the Green’s function 
for the continuous-state random walk and approximate it by its continuous counterpart. 
In Section [3] we give finer estimate for the Green’s function and prove Proposition [l] In 
Section H] we prove (JH). 

In the Appendix, we prove the asymptotics for the potential kernel of the continuous- 
state random walk. 


2. Preliminaries 


(9) 


2.1. Continuous-state random walk. Recall that our continuous-state random walk 
is defined by S n — S 0 + X 3 + X 2 + ■ ■ ■ + X n for n G N U {0}, where AYs are i.i.d and 
each Xi is uniformly distributed in the disk of radius h > 0. We use P x to denote the 
conditional distribution of {S n } n > 0 , given S 0 = x\ and write E x for the corresponding 
expectation. For simplicity we suppress the h dependence in P x and E x since all random 
walks have h dependence in this paper. We first define the transition density for such 
random walks: 

p(0,x,y) = S(y-x), 

P(k,x,y) = lim 40 pX({S *~ 2 yl - e \ k > 1, 
where 5(y — x) for fixed x is the delta function giving unit mass to the point x. For a 
simply connected and bounded domain D , we define the first time the continuous-state 
random walk leaves D as To, i.e., 

T d := inf {k > 0 : S k <£D}. 

We denote the transition density for continuous-state random walk killed on exiting D 
by pd, i-e., for x G D and y G C 

p D (0,x,y) = 8(y - x) 

p D (k, x, y) = lim 40 k<r D ) ^ k y 1 

One immediate consequence of the Markov property for S n is for y G D 

[ PD(k - l,x,y + £)d£, k > 1. (11) 

<B(0,h) 


( 10 ) 


p D (k,x,y) = 


nil 2 


Note that the same equality holds if pr> is replaced by p in (ITT]) . The discrete-time Green’s 
function for D is defined by 

OO 

G h (x,y) = ^2p D (k,x,y), x G D, y G C \ {a:}. 

k =0 

Recall the definitions of Di, D 2 and ZR in the introduction. 

Lemma 1. For a fixed x G D, the discrete-time Green’s function satisfies 

A h G h (x, y) = 0, y e D\{x} 

G h (x,y) = 0, y e D 3 . 


(12) 
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Proof. The boundary condition is obvious. A h Gh(x,y) = 0 for y G D \ {x} follows from 
(TTTD and the Fubini-Tonelli theorem. □ 


The following is a maximum principle for the discrete Laplacian defined by (171) . 

Lemma 2. Let D be a simply connected and bounded Jordan domain. Suppose for some 
M > 0, f is a function in D U D 3 satisfying 

\A h f(z)\<M,zeD. 


Then we have 

sup |/(*)| < 2Mr 2 (D)/h 2 + sup |/(*)|, 

z£D ZGD3 

where r(D) := inf{r > 0 : DUD^ C 5(0, r)} is the radius of the smallest circle (centered 
at 0) circumscribed about hUh 3 . 

Proof. The proof is similar to the proofs for Lemmas 23.4 and 23.5 in [5|. □ 

Next, we write the discretization error fh — f in terms of the discrete-time Green’s 
function. 

Lemma 3. Under the assumption of Proposition^ we have 

A(0) - /(0) = f G h (0,z)A h f(z)dz. (13) 

Jd 2 

Proof. We use a martingale formula (see Lemma 1 of [9]). We extend / to D U D 5 by 
setting f(P) = f(P) for P G D$ where P = arg rniri{|z — P\ : z G dD}. It is clear that 


TO—1 

f(S n )-J2 A hf(S k ),n>0 

k =0 

is a martingale with respect to T n := cr(5o, S\, ■ ■ ■ , S n ) (note that this is actually true for 
any /). Then the optional sampling theorem and the dominated convergence theorem 
(noting that / is bounded) give 


To-1 

B[/(StJ-/(S„)] = B[^ A„/(S t )]. 

k =0 

Suppose = 0. Then we have 


A(o) - /(o) 


T]j — 1 oo 

E°{J2 A, f(S t )] = £»[ J2 A h f(S t )I {tiTD - 1} ] 

k =0 k=0 

oo « 

V / p D (k,0,z)A h f(z)dz 

k =o Jd 


G h (0,z)A h f(z)dz 


id 



G h (0,z)A h f(z)dz 


where we have used Fubini’s theorem (since E°\ J2k=o 1 ^/(SV-)! A 2||/|| oo 5 0 T£ ) < oo) 
and the fact that A hf(z) = 0 for z G D\ (by the mean value property for /). □ 
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In order to estimate the discrete-time Green’s function, we need some result concerning 
the potential kernel. For igC, let a,h(x) be the potential kernel for the continuous-state 
random walk defined by 

OO 

Uh(x) = (14) 

k =1 

Note that we do not include the k = 0 term in the above sum so that a,h(x) is a function 
rather than a distribution. For simplicity, we write a(x) for ai(x). So a(x) does not 
depend on h. 

Lemma 4. The series in m is convergent, i.e., a(x) is well-defined. A h.a(x/h) = 
4/ B(Q,h){x), and there exists a constant Co > 0 such that 

4 

a(x) = — In |x| + C 0 + 0(|:r|~ 2 ) as x —> oo. 

7T 

Since the analogous result is standard for discrete-state random walks but the proof is 
quite long, we will present the proof in the Appendix. 

Remark 5. See [6] for a proof of the above lemma for discrete-state random walks. See 
also [II] for a different proof. [3] has a proof for a discrete-state random walk on isoradial 
graphs. 

2.2. The difference between discrete-time and continuous Green’s functions. 

For the domain D defined in Theorem [fl the continuous Green’s function Gd(0, z) is 
the unique harmonic function on D \ {0} such that Gd(0,z) —» 0 as z —> dD and 
Gd( 0, z) = — In \z\/(2it) + 0(1) as z —» 0. Then Gd(0, z) — — In |z|/(27t) + if(z) where 
if(z) is defined by ((171) below. Following the estimate of discrete Green’s function for 
simple random walk in [15j, we prove the following estimate for the continuous-state 
random walk. 

Lemma 5. Suppose D C C is a simply connected and bounded Jordan domain, and dD 
is analytic. Assume 0 G D. Then there exists a C > 0 independent of h such that 

\h 2 G H (0,z)-8G D (0,z)\<Ch (15) 

uniformly for z E D satisfying \z\ > \[h. 

Remark 6. See Theorem 1.2 of [10] and Appendix A of [Ij for results on the difference 
between the discrete Green’s function for the simple random walk and the continuous 
Green’s function for domains without any smoothness assumption on the boundary. 

Proof. We apply the method introduced in [T5]. The basic idea of the proof is that we 
will find the relationship between Gh(0,z) and Gd(0,z) via the potential kernel a(x). 

Without loss of generality, we may assume B(0,Vh) C D. For z E C, let H h {z ) = 
—h 2 5(z) — a(z/h ) — 81nh/(27r) + C 0 . Then by Lemma [H, we have 

H h (z) = —8 In |z|/(27t) + 0(h 2 /\z\ 2 ) for z E D \ B( 0, Vh). 

For z G C\{0}, let Ch(z) := h 2 Gh( 0, z) — H h (z). Note that both terms in the difference 
contain a delta function at 0, but these delta functions cancel. We have 

OO 

6h(z) = h 2 'y^fpnjk, 0, z) + a(z/h ) + 8 In h/ (27r) — Cq, z G C \ {0}. 

k =1 
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We can use the equation above to define e^(0). Then by (112]) . Lemma HJ and the above 
discussion, we have 


/ A h e h (z) = 0, z G D 

\ e h {z) = 8In \z\/( 2 ti) + 0(h 2 /\z | 2 ), z G D 3 . 

Suppose 'if(z) is the harmonic function of z G D satisfying 

f A ip(z) = 0, z G D 

\ if(z) — In |z|/(27t), zeOD. 


(16) 


(17) 


Note that 'if(z) = Gd(0,z) + ln|z|/(27r) for z E D \ {0}, and Gd(0,z) can be extended 
to a harmonic function in a domain containing D \ {0} (see Lemma II.2.4 of [7]) and 
In \z\/(2n) is harmonic in a domain containing D \ {0}. Hence if(z ) can be extended to 
a harmonic function in a domain containing D U D 3 when h is small. Therefore for all 
small h > 0 

f A h i/j(z) = 0, zed . . 

\ 'ifj(z) = In\z\/(2n) + O(h), z G D 3 

where the O(h) comes from the harmonic extension of ip(z). Subtracting 8xfrom 
(Hi, we get (note that we assumed B(0, Vh) C D ) 


f Afc[e k (z)-8^(*)]=0 z G D 
\ e h (z) - 8^(z) = O(h), z e D 3 . 

Then the maximum principle for Ah, Lemma El implies 


(19) 


eh(z) — 8 i)j(z) = O(h), z e D. 


Therefore, 

h 2 G h (0, z) = e h (z) + H h (z) — —8 In \z\/ (27t) + 8'i/j(z) + 0(h) 

= 8G/)(0, z) + O(h) 

where the second equality is true if \z\ > Vh. 

□ 


3. Proof of Proposition dj 


3.1. Further estimate of the discrete-time Green’s function. Recall that if the 
boundary of D is analytic, then the Poisson kernel H d (0,x) for 0 G D,x G dD is defined 
by 


H d (0,x) = 


dG D (0,x) 


On, 


where Gp is the continuous Green’s function and n^, 
have the following easy estimate 


is the inward unit normal at x. We 


Lemma 6. Suppose D C C is a simply connected and bounded Jordan domain, and OD 
is analytic. Assume 0 G D. Then for l e [0, h\ 

G d ( 0, x + ln x ) = IH d ( 0, x) + 0(h 2 ), 

uniformly for x G OD . 

Proof. Note that Gd( 0,x) = 0 for any x G OD. Lemma II.2.4 of [7j implies Go could be 
extended to a harmonic function in a domain containing D \ {0}. The lemma follows by 
Taylor expansion of Gd(0, x + ln x ) about x with coordinate directions n x and — in x . □ 
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Next, we improve the estimate in Lemma 0 We need a Beurling-type estimate for the 
continuous-state random walk. For u G D and V C dD, we define 

dist^i (u; V) = infji? > 0 : u and V are connected in D fl B(u, R)}. 


Lemma 7. There exist two absolute constants (3 > 0 and C > 0 such that for any simple 
connected domain D , u G D and V C dD we have 


P u (S Td eV)<C 


dist(w, dD) 
dist^ (w; V) 


P 


1 


where St d is the orthogonal projection of St d onto dD, i.e., St d = argmin{|z — St d 
z G dD}. 


Proof. The proof is almost the same as the proof of Proposition 2.11 in [[3j for a discrete- 
state random walk. □ 


The following proposition is an improvement of Lemma [5] for points very close to the 
boundary of the domain. 

Proposition 2. Suppose D C C is a simply connected and bounded Jordan domain, and 
dD is analytic. Assume 0 G D. Then for any e G (0,1/2), x G dD and l G [0, h\ 

h 2 G h (0, z) - 8 G D ( 0 , 2 ) - 8H D (0,x)E li [\lm(S Tn )\] = 0(h 1+ ^) + 0(h 2 ~ 2e ) 

where z = x + ln x G D 2 , and the big O terms only depend on D. 

Proof. From the proof of Lemma 0 we know it is enough to prove for z = x + ln x G D 2 

e h (z) - 8if(z) = 8H D (0,x)E H [\lm(S TB )\] + 0(h 1+ ^) + 0(h 2 ~ 2e ). (20) 

If we write out the 0(h) in (Hi?]) , we see eh(z) — 8 'if(z) satishes 

( A h [e h (z) - 8 ^( 2 )] =0 z E D 

\ e h (z) - 8ip(z) = —8G d (0, z) + 0(h 2 ), z G D 3 , 

where we have used the same notation (i.e., Go) for the harmonic extension of the 
continuous Green’s function. It is easy to see that Lemma [ 6 ] is also true for l G [— h, 0], 
i.e., x + ln x G D 3 . Let Fh(0,z) be the solution of the following discrete Dirichlet problem 

( A h [F h (0,z)] =0 z G D 

\ F h ( 0, z) = -IH d ( 0, x), z = x + ln x G D 3 , 

Then (1201) follows from Lemmas [2] and [Gj and the following claim: for 2 = x + / n ; „ G D 2 

F h (0,z) = H D (0,x)E l %lm(S TB )\] + O(h 1+ ^) + O(h 2 - 2e ). (21) 

Let t x be the unit tangent vector at x G dD. One can choose either = in x or t x = —in x . 
Let 7 (s), 0 < s < \dD\ be an arc-length parametrization of dD. For any x G dD, let 
a(x) G [0, |<9-D|] such that 7 ( 0 / 2 ;)) = x, i.e., a is the inverse of 7 . Since dD is analytic, 

d// ( 7 [c r (ic) — h}~ e ,a(x) + • [—h l ^ e ,h 1 ~ e ]) < Mh 2 ~ 2e 

uniformly for x G dD, where dn is the Hausdorff distance and t x ■ [—h 1 ~ e ,h 1 ~ e ] is the 
tangent line segment centered at x with total length 2 h v ~ e . For any z 3 = X\ + /in xi G D 3 
where 27 G dD, dehne x : D 3 —> dD and l : D 3 -E [— h, 0] such that 2 /zi) = 27 and 
l(z\) = l\. Note that x(z\) and l(z\) are uniquely dehned if h is small enough. For any 
Zq = Xq + Iq 1 lI Xo G D 2 , the comment after (JH]) implies 

^( 0 , 20 ) = E Zo [-1(S Td )H d (0,x(S T d ))]. 








Lemma d gives (noting that x(St d ) = St d ) 

P Xo+lon *o (x(S T d ) i l[a{x Q ) - h}-%a(x o) + h 1 - 6 ]) = 0 (/^). 

Therefore, 

-Pft,(0,^o) = -E'" 0 [—■ 1(St d )Hd( 0, x(ST D ))I{ 1 [a(x 0 )-h 1 -^cT(xo)+h 1 -^}}(x(ST D ))} +0(h 1+€/3 ). 

The smoothness of Ho( 0,x) in x G dD implies Hd(0,x(St d )) = Hd(0,xo) +0[h 1 ~ e ) 
if x(St d ) G 7 [cr(o;o) — /r 1_e , cr(xo) + h l ~ e ]. Hence 

Fh(0,z o ) = i/ I) (0,xo)^ 2o [-/(5T D )/ { ^(xo)-^-s ff (xo) + ^}(^( 5 '^))] +0(/r 2 - e ) + 0(h 1+ ^). 
Let 

B X0 {h}~\ Mh 2 ~ 2e ) := t X0 • [-h 1 ” 6 - Mh 2 ~ 2e , h 1_e + Mh 2 ~ 2e ] + n Xo • [-M/i 2_2e , Mh 2 " 2e ] 
be the rectangle centered at Xq with length 2 h l ~ e + 2 Mh 2 ~ 2e and width 2 Mh 2 ~ 2e . Then 
7 [cr(a;o) - h^^afx 0 ) + /i 1-e ] C B Xo (h}~ e ,Mh 2 ~ 2e ). 

It is not hard to see 

P Z0 (S(T d ) G B^-^Mh 2 - 26 )) < sup P [i (|Im(S' TH )| < 2 Mh 2 ~ 2e ) = Oih 1 -^). 

0<l<h 

Let P^ o (h 1_e , Mh 2 ~ 2e ) := C \ B XQ [h}~ e , Mh?~ 2e ). Then we have 

Ffi(0, ^o) Bf)(0,X())E [ )I{'y[<T(xo)—h i -~ e ,<r(a:o)+/il— e ]} ) ) 

x/B go ( fc i-. w - 2 «)(5(T ZJ ))] + 0 (h 2 ' 2e ) + 0 (/i 1+ ^) 

= P D (0,x 0 )P a [|Im(P TH )|] + 0(h 2 ~ 2e ) +0(h 1+e P). 

This completes the claim (and hence the proposition) since z 0 is arbitrary. 

□ 


3.2. A change of variables formula and an estimate of A^f. From Lemma El we 
know the difference of //,(()) and /(0) can be represented by a two-dimensional integral 
in D 2 . We give a change of variables formula for such an integral in the following lemma. 

Lemma 8. Suppose D is a simply connected and bounded Jordan domain, and dD is 
analytic. Let F be a Lebesgue measurable and bounded function on D 2 . Suppose z(t ) = 
(u(t),v(t)), 0 < t < 1 is a parametrization of dD such that u' 2 (t) + v' 2 (t) ^ 0 for any t. 
Then the following holds for all small h > 0 

J D F ( z ) dz = J o J o F ( z ( t ) + in *( t ))( 1 - 1 ^ u ,2 + v , 2 y / 2 ) dl ^ u ' 2 F v ' 2dt ^ 

where n z (t) is the inward unit normal at z(t). In particular, for all small h > 0 
[ F{z)dz = (l + 0(h)) f [ h F{£ + lnt)dl\dt\, 

J D 2 J dD JO 

where 0(h) only depends on D. 


Proof. The lemma follows from the change of variables 

lv' 


z = (x,y) = z(t) + ln z( t) = ( u(t ) 


,,v(t) + 


lv! 


r)- 


yju ' 2 + v ' 2 ’ yju ' 2 + v ' 2 ‘ 

Note that one needs to pick h small to make sure the above change of variables is one- 
to-one. □ 


We still need some estimate on A hf for the / defined in Proposition [H 
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Lemma 9. Under the assumption of Proposition d for any x G dD and l G [0, h] we 
have 


1 r)f(r) 2 _ l 2 

A h f(x + ln x ) = —— ' 2 [~/? 2 \/h 2 - l 2 + — \//i 2 - l 2 - lh 2 arccos(Z/h)] + 0(/z 2 ), 

nlr on x o o 

where 0(h 2 ) only depends on f and D. 

Proof. 


A h f(x + ln x ) = —— / [/(a; + Zn x + 0 - f(x + ln x )}df 

•2n 

r[f(x + Zn x + re* 0 ) — f(x + Zn x )](Z#(Zr 

r[f(x + Zn x + re* 0 ) — f(x + ln x )]d0dr (22) 

Since / is harmonic in D, the first integral in (j22j) is zero due to the mean value property 
for /. Since dD is analytic, there exist a conformal map T from D to D and an e > 0 
such that T can be extended to a conformal map of (1 + e)B. This implies that the arc 
dD D B(x, h ) can be approximated by the tangent line segment at x with fixed length of 
0(h). More precisely, there exist constants C\ > 0 and C 2 > 0 such that 

dD fl B(x, h ) C t x • [— Cih, C\h] + n x ■ [—C 2 h 2 , C 2 h 2 ] 




uniformly for all x G dD, where t x is the unit tangent vector at x. Therefore (HI2D gives 


A h f(x + ln x ) = 
1 


| nh /»7r/2+arcsin (l/r) 


nil 2 

h /»arccos (l/r) 



r[f(x + ln x + rn x e* 0 ) — f(x + ln x )]d6dr 


l — 7r/2—arcsin(Z/r) 


7T ll 2 



r[f(x + ln x — r n T e* 0 ) — f(x + ln x )]d6dr + 0(h 6 ) 


'l J — arccos(Z /r) 

•— I\(%■) l ) + I 2 (#> l) + 0{h?'). 


Noticing that / G C 2 (D), by Taylor expansion of / about x with coordinate directions 
n x and — m x , we get 


1 /-h />7r/2+arcsin (l/r) f) f ( r ) 

h(x,l ) = —-J / / r[——-(rcos6>)]cZ6><ir + 0(h 2 ) 


Ji 7r/2—arcsin(Z/r) 


where we used the fact the coefficient for (which is rsin#) is an odd function of 6. 

Similarly, by the definition of / in D 3 and the Taylor expansion, we have 

•1 ph /»arccos (l/r) 

h(x,l) = —rW / r[f(x-in x rsm(9))-f(x + ln x )]d9dr + 0(h 2 ) 

Ji arccos(//r) 

1 fh /-arccos (Z/r) flff '* 

1 1 r[-l^l}d6dr + 0(h 2 ). 


nh 2 



l J — arccos (Z/V) 


<9n. 


The lemma follows by simple computations of Ii(x,l) and I 2 (x,l). 


□ 
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3.3. Proof of Proposition [Q Now we are ready to prove Proposition [0 

Proof of Proposition d Recall that J g(£,)oJh( 0, |d£|; D) = fh{ 0) and f g(f)co( 0, \df\] D ) = 
/(0). By Lemmas [3] and El we have 


= a + 0(h)) 

Using Lemma [9J we get 


A(0 )-/(0)= / G h (0,z)A h f(z)dz 
Jd 2 

/ G h (0,x+ln x )A h f(x + ln x )dl\dx\. 


'8D JO 


D{h) : = 


'8D JO 

1 f 


Gh( 0, x + ln x )A h f(x + ln x )dl\dx\ 
df(x) ' h 


7r 


h 2 


dn 


G h ( 0,a; + ln x ) 


dD ull x Jo 


2 _ 7 2 

*[-/r\// 7, 2 — l 2 H- \/h 2 — l 2 — III 2 arccos(//h) + 0 (h 4 )]d/|dx|. 

3 3 


Applying Proposition [2j we obtain 

s /' df(x) f h 


D(h ) = 


dD Jo 


(G d (0 ,x + InJ + He>(0, x)£ ,/ *[|Im(S , T II )|] + o(h)) 


7r/i 4 ^ 

*[^-h 2 Vh 2 — l 2 + l—Vh 2 — l 2 — Ih 2 arccos(Z//r) + 0(/i 4 )]<i/|ch| 
o o 


Substituting the estimate in Lemma El into the above equality, we see 

s r df(x) ' h 


D(h) = 


nhiJ g ni J (lH D (0,x) + H D (0,x)E h {\lm(S Ta )\] + o(h)) 

2 _ /2 

[—/? 2 \/ h 2 — l 2 H-a//i 2 — l 2 — Ih 2 arccos(Z/h) + 0 (Z 7 4 )]cZZ|(i:r|. 

3 3 


By the change of variables Z = h cos 9, we have 


16 8 


t/2 


Z?(/i) = + A / (sin 2 ^ — (sin 4 0)/3 — 6*cos6 l sin6 l )U^“ 1 se (|Im(S'T II )|)d0] 

457T n l n 


* [ df(x) 

J8d dn x 

which completes the proof. 


H D ( 0, x)|dcc| + o(h), 


□ 


4. Proof of density for limiting measure 


Proposition 3. Under the assumption of Proposition^ let if be a conformal map from 
D to D which sends 0 to the origin. Ford e [0,27 t], define m(9) = |^ / (^ _ 1 (e* e ))|. Then 
we have 



H D (a,z)UA\ dz \ 

on. 



(goif l )(e^)p((f)dcf , 


where 


p(<P) 


1 f 2n m(9) — m(<f) — m!((f) sin( 6 * — <f) 

4n 2 J 0 1 — 005(6* — ( f ) 


li 



















Proof. We define 


I (g)=[ H D (0,z)-jr^-\dz\. 

JdD 9 n, 

Using the change of variables 0(z) = e ld for z e dD, we have 

d f( z ) = ° *P~ 1 )(e ie ) _ , Q , Um (/ ° ^~ 1 ){e i6 + en el9 ) - (/ o if- 1 )^ 6 ) 

dn z dn pi e 40 e 


= 777 (d) lim 
40 


(/°0 x ) ((l-e)e i0 ) - (foil) l )(e ie ) 


Under the same change of variables, the harmonic measure Ho(0, z)\dz\ transforms to 


±dd. So 


-2 tt 


1(g) = — / m(9) lim 


(/ 0-0 *) ((l-e)e i0 ) - (/°0 l )(e ie ) 


2vr 


40 


dO. 


By the assumption of the proposition, m(6) is a smooth function of 9 and /o-0 1 e U 2 (D). 
So the mean value theorem and the bounded convergence theorem give 


»2tt 


I (9) = 7T lh ff 

Z7T e4,0 J q 


(/ O -0 x ) ((1 - e)e ie ) - (/ o 0 ^(e* 0 ) 


dO. 


Note that / o -0 1 is the harmonic function on D with boundary data g o 0 1 . We now 
have 


r-2-n- 


7 (. 9 ) = hu lir 0 / m ( 0 ) 

Z7T e.J.0 /n 


((1 -e)e i 0 ,e^) [( 90 ^ ^(e^) - (# o 0 ^(e* 0 )]^ 


dd. 


Fubini’s theorem and interchange of 6 and 0 imply 


-2 IT 


f27T 


/ (^) = 1 4o/ 0 (3°^ ^(e^j 

Reflection symmetry implies 


m(6 l )ffB ((1 — e)e 20 , e l<?i ) — m((j))H® ((1 — e)e*^, e 10 ) 


27T6 


dOdcj) 


r 2 -k 


770(0) sin(d — 0 )/Jb ((1 — e)e u e* 0 ) dd = 0. 


So we can rewrite our integral as 


r-27T 


1(a) = —lim / 
27T 40 J o 


(s »r‘)y) i-4“,e‘*) 


m(0)i/B ((1 — e)e®^, e* 0 ) m'(0) sin(d — 0)dd D ((1 — <f)e l<t> , e* 0 ) 


]ddd0. 


Rotation and reflection symmetries imply 

tfB((l-e)e i0 ,e^) = H n ((1 - e)e^, e i0 ) = H 0 ((1 - e), e^) 


2e — e 2 


27T 2 — 2 008 ( 6 * — 0 ) — [2 — 2 cos(d — 0 )]e + e 2 ’ 

where we have used the explicit expression for the Poisson kernel in D (see, e.g., Example 
2.16 of [ 12 ]). So we now have 


1 


r2n 


-2 IT 


/(9) = 3 v^L (9 ° ,/rl)(e ' : * ) 


(2 — e) [ 777 (d) — 777(0) — 777'(0) 8111(6* — 0)] 

2 — 2 cos(6* — 0) — [2 — 2 cos(6* — 0)]e + e 2 
12 


dOd<f>. 

















Viewing 2 — 2 cos(0 — 0) — [2 — 2 cos(0 — 0)]e + e 2 as a quadratic function of e, it is easy 
to show for any 0 < e < 1 one has 

2 — 2 008 ( 6 * — (f>) — [2 — 2 cos(6* — 0)]e + e 2 > [1 — cos(0 — 0)]/2. 


Therefore 

(2 — e)[m(0) — m(0) — m'(0) sin(0 — 0)] 4|m(0) — m(0) — m'(0) sin(0 — 0)| 

2 — 2 cos(0 — 0) — [2 — 2 cos(0 — 0)]e + e 2 — 1 — cos(0 — 0) 

L’Hopital’s rule applied to the left hand of side fl23|i implies that it is a bounded function 
of 0 G M and 0 6 I (using the periodicity). So by the bounded convergence theorem 
(noting that g o -0 -1 is also bounded) we have 


r-2-K 


1(9) = 


4tt 2 

which is the desired result. 


(c/°0 )(e 



m(0) — m'(0) sin(0 — 0) 
1 — cos(0 — 0) 


dddcj), 


□ 


Let us remark that Theorem [Tj follows from Propositions Q] and [3l and the change of 
variables z = 0 _1 (e*^). 
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5. Appendix 

In this appendix, we prove the asymptotics for the potential kernel described in Lemma 
EO We follow the methods introduced in Section 12 of ca and [6]. Let 0(0) be the 
characteristic function of the continuous-state random walk with h = 1, i.e., 

0(0) = Ee iX ' e 

where 6 = (0i,02) and X = (AW, X^) is uniformly distributed in the disk of radius 1. 

Lemma 10. 

m = i-^+)§)+o(i9iv^o 

| <p(0)\ < 1 min^iC 1 , \0 2 \-'},0 -f oo. 

IT 

Proof. The hrst estimate in the lemma follows by Taylor expansion, while the second 
follows since 

1 r 1 /V 1 -*? 4 

|0(0)| = |— / / _cos(xi0i) cos(x29 2 )dxidx2\ < — r^ 

71 J-i J-yf i-xl 7r |021 

and the symmetry of 0 1 and 02 - □ 

The following lemma says our potential kernel is well-defined. 

Lemma 11. 

Tl 

a(x ) = lim y \p(k, 0, 0) — p(k, 0, x)] 

71—>-00 ' J 

k=1 

2 1 /* 1 iO-x 

= 5Z[p(fc,°,°) ~p{k,0,x)} + / - e — --0 3 (0)d0 

( 2?r ) 0R2 1 - 
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Proof. By applying the continuous inversion formula, the proof is similar to the proof of 
PI in Section 12 of [14] if one can show 0 3 (d) £ P^M 2 ). The latter is true because 
of Lemma US □ 

Let Q{9) = E{X-Of 2 = and 0(d) = 1/(1 —0(d)) —2 /Q{9). Then Lemma ITOl implies 

-0(d) = 1/3 + 0(|d| 2 ) as d —> 0; |0(d)| < 2 as d —>■ oo. (24) 

Now we have all ingredients to prove Lemma |4j 


Proof of Lemma [/} By Lemma Qj] and the evenness of 0, we see that 

a ( x ) = ^[P( k ,0,0)-p(k,0,x)} + \ I -— g ^ 0 3 (d)dd 
k= i ^ R2 I I 

+ 77/12 / t 1 - (25) 

(2tt) 2 y R 2 

By the estimates in Lemma dS and (HTH) . and the Riemann-Lebesgue lemma, 


70- [ (1 - e »*)4ie)<fH0)de -> —!_ [ 4,(e)<j?(e)de as | x | -> 00 
(2vr) 2 ,/ K 2 (2 tt) 2 V 

which is a constant contributing to Co in the Lemma. 

This gives the first o(l) term 

- / e* e ' 3: 0(d)0 3 (d)dd = -—^-7 [ cos ( x ■ d)0(d)0 3 (d)dd. (26) 

(2vr)- J R 2 (27t) 2 J R 2 

Let 5 := 5(0,7r) := {z : \z\ < 7r} and B c = M 2 \ B. Then the first integral together 
with the attached multiplicative term in (j25|) can be written as the sum of the following 
two integrals 


IAx) ■= 


h{x) := 


7 


7T 2 / 

" Jb c 


1 — cos(a: 

I d | 2 

'^0 3 (d)dd 

(27) 

1 — cos(a; • 
|d| 2 

% 3 (d)dd. 

(28) 


By the estimate in Lemma US and the Riemann-Lebesgue lemma we have 

T / \ 2 /• 0 3 (d) 

H x ) -»• ^ / “mW" asx^oo, 


7T 2 / 

" JBC 


|d/ 


which leaves the second o(l) term 


cos(a; • d) 3 


tt 2 J Bc |d 2 


0 3 (d)dd. 


We rewrite h(x) as follows 

AW = 4 f lz^l de+ l f * , W - 1 <w+ . 2 


7T 


/ B 


Id/ 


7T 


IB 


Id/ 


7T 


IB 


cos(x • d) 


(29) 


(l-0 3 (d))dd. (30) 


By Lemma [TS the second integral in (13UD is a constant contributing to Co in the lemma, 
and by the Riemann-Lebesgue lemma the last integral in (1HUD gives the third o(l) term 


7 r z 


cos (a; • d) 


(1 - 0 3 (d))dd, 
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(31) 

















and the first integral together with the attached multiplicative term in (150]) is equal to 
(using the proof of P3 in Section 12 of [H]) 



(32) 


where 7 is the Euler’s constant. 

It is clear that 7 + lm/r + In \x\ + ln(sino;) in (152]) as a function of a is integrable from 
0 to 7t/2, so the fourth o(l) term is 


r w/2 


cosn 2 

- duda = — 


7T 


'7r|a:| sina 


U 


7T- 


' B c 


cos(x • 6) 


dO. 


(33) 


where the equality follows by reversing the procedure which led to ([32]) . 
Adding the four o(l) terms, i.e., d26]) + (l22]) + (l3l]) + (l33]) . we get 


(2*0 


cos(x • 9)'4t{9)(f^{9)d9 H—- 


7T 


cos(x • 6) 


(1 ~(t?(6))de 


4 n 2 


cos(x ■ 6) [—— — 


\6) 


101 2 1 -m 


}d9. 


(34) 


Noting that cos(x • 9) = V • b(0) where b(0) = sin(x • 9){x i/|x| 2 , X 2 /|x| 2 ), the divergence 
theorem gives 


cos (z-0)[j^ 

. 8 

W 


\9) 


lim / 

N ^-°o JB(0,N) 

lim lx A 1 

N —>00 


1 - 4(9) 

1 


]d9 

V ■ b {9)dJ9 


*B(0,N) 


X2_. . 8 (j) 3 (9 ) 

’|x|’|x| J Vl 2 1-0(0) 


sin(x ■ 9)d9 


(35) 


We can apply the divergence theorem again to (1351) . As a result, we see that (134]) has 
order 0 (|x|~ 2 ). 

Therefore, the proof of Lemma 0] is complete if one can show A h.a{x/h) = y/_B(o,/i)(x). 
But the latter is easy to verify (note that a(x/h) = h 2 dh(x)). 

□ 
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